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In this chapter we discuss some possible physical pictures that de- 
scribe the constitution of the inner crust of compact objects. Different 
relativistic models both with constant couplings and density dependent 
ones are used. We calculate the liquid-gas phase transition in asymmetric 
nuclear matter from the thermodynamic and dynamic instabilities. The 
C^ ■ equations of state used to describe the crust are related to the crust-core 

transition properties. Cold and warm pasta phases with and without al- 
pha particles are constructed. The influence of the pasta phase and its 
internal structure on the diffusion coefficients associated with Boltzman 
transport equations used to simulate the evolution of protoneutron stars 
are shown. Finally, the possible existence of bare quark stars and the 
effects of strong magnetic fields on quark matter are considered. Open 
questions are pointed out. 
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1. Introduction 

The internal constituents of compact stars are a great source of speculation ifTTl . 
They could be made of hadronic matter only, quark matter only (a possibility 
arisen by the Bodmer-Witten conjecture [2,j3J), or they could be hybrid. Hybrid 
stars may have in their interior hadrons and quarks with or without a condensate. 
If the star is composed of quark matter only, it may be a bare star and we will 
tackle this point at the end of this chapter. On the other hand, if it is made of 
hadrons only or if it is a hybrid object, it has a crust, our main interest in what 
follows. The crust mass and thickness depend on the equation of state (EOS) 
that describes the star. They also depend on the total stellar mass; an increase 
of the stellar mass increases the gravitational pull within the crust resulting in 
the thinning of the crust. An estimation of the size of the crust done in [ 5 1 gives 
values varying from 1.01 to 0.29 Km. In particular, we will discuss the possible 
existence of inhomogeneous structures, i. e. the pasta phase, in the inner crust 
due to the competition between the strong and Coulomb interactions and their 
implications for the properties of the crust. 

Through out this chapter, we consider a system of protons and neutrons 
with mass M interacting with and through an isoscalar-scalar field cj) with mass 
rus, an isoscalar- vector field V^^ with mass m.^, an isovector- vector field b^ 
with mass nip. We also include a system of leptons composed by electrons and 
muons, electrons and neutrinos or just electrons, depending on the problem. The 
Lagrangian density reads: 

i=p,n l=e,u,fi 

where the nucleon Lagrangian density has the form 

Ci = tl;i[-f^iD^-M*]tl;i, (2) 

with 

iD^' = i^^'-ryV>'-'^T■h^'-e^-^A'', (3) 

M* = M-Tscl). (4) 
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The lepton Lagrangian density is given by 

A = V^z[7/.(^9^ + e^'^)-mz]^z, (5) 

where e, mi stand for the charge and mass of the lepton, respectively, and the 
meson Lagrangian densities are 



1 / 1 



2 V 2 



%M'"' + mlV^V 



l^P = U-^B.^-B^'^ + mlh.-h^ 



C — —tp p^^>^ 



where 9.^^ = d^V^ - d^V^ , B^j, = d^h^ - d^h^, - Tp{h^ x b^) and 
Fiiv = dpAy — duA^. The parameters of the models are the nucleon mass 
M = 939 MeV, the coupling parameters Tg, T^, Vp of the mesons to the nucle- 
ons, the electi^on mass rrif. = 0.511 MeV, the muon mass nip = 105.66 MeV 
and the electromagnetic coupling constant e = y^47r/137. The electron neu- 
trino mass is considered to be zero. In the above Lagrangian density r are the 
Pauli matrices. When density dependent models are used, the non-linear terms 
are not present and, hence, k = A = and the density dependent parameters 
are chosen as in ||4l|6l|7l. For the parametrizations with constant couplings, Fj 
is replaced by gi, where i = s, v, p as in the NL3 parameter set HI. 

Some expressions are often used next. Yp and Yl refer to the proton and 
lepton fractions respectively. If matter in /3-equilibrium is considered in a sys- 
tem of protons, neutrons, electrons and possibly trapped electron neutrinos, one 
has: 

Hp = Hn- ilJ'e- fJ-u)- (6) 

For neutrino free matter fi^ = 0. Neutrality of charge requires pp = pe + Pp 
and, if pp is present, its chemical potential is equal to the electron chemical 
potential. 

2. Instabilities in Nuclear Matter 

It was shown in |l9l[l3 that the pressure and density at the inner boundary of the 
crust (transition pressure and transition density) define the mass and moment 
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of inertia of the crust. This establishes a relation between the equation of state 
(EOS) and compact-star observables. In this section we show how an estimation 
of the transition properties may be obtained from the thermodynamical binodal 
and spinodal surfaces, or the dynamical spinodal. A better estimation, obtained 
from the pasta phase calculation, will be discussed in section [3l 

The liquid-gas phase transition at subsaturation densities is a well known 
feature of the nuclear EOS [11] . It corresponds to the presence of a negative 
curvature of the free energy density T . In this case the system is unstable against 
separation into two infinite homogeneous phases. The spinodal surface, which 
limits the unstable region in the (T, pp, pn) space, is defined by the cancellation 
of the determinant of the free energy curvature matrix [.12.1 : 

""U^/^Jt' lis i^j- 

The eigenvalues of the curvature matrix are given by 

A± = ^ (Tr(C) ± VTr(C)2 - 4Det(C) 

and are associated to the eigenvectors 6p±. Inside the spinodal region the low- 
est eigenvalue A_ is negative. The direction of the vector Sp- defines the direc- 
tion of instability, which generally dictates a distillation effect corresponding to 
a phase separation into a high density symmetric matter and low density neutron 
rich matter [121 . 

In the left panel of Fig {T] the spinodal sections obtained for different rela- 
tivistic mean field (RMF) models are plotted in terms of the total density p and 
the asymmetry parameter 

6 = iPn- Pp)/P = 1 - 2Yp. (7) 

While for symmetric matter most of the models show a similar behavior, at high 
asymmetry and/or temperature models differ. For this reason, establishing con- 
straints on the equations of state based on experimental results is an important 
step, and a discussion on this aspect is done in section |4] In the right panel of 
Fig. [H the spinodal sections for several temperatures and two models (NL3 HI 
and TW [4J) are an example of the possible existing differences: for TW with 
a symmetry energy and corresponding slope equal, respectively, to 32.76 and 
55.30 MeV, the transition density (identified by a full dot) occurs at a larger 
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density and for a smaller proton fraction than for NL3 with a symmetry en- 
ergy and slope equal, respectively, to 37.34 and 118.30 MeV. Also, for TW the 
spinodal surface extends to a larger temperature. 
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Figure 1 . Left panel: the spinodal section for different RMF models models in 
the p — S plane. Fig. taken from (T\. Right panel: the spinodal section for NL3 
and TW and different temperatures. The /3-equilibrium EOS at T = (T = 10) 
MeV is represented by a full (dashed) line. The full dot identifies the transition 
crust-core at T = 0. Fig. taken from llT2l . 



Instead of using a thermodynamic approach, in 11131 the transition density 
was estimated in a local equilibrium approximation by calculating the density 
for which matter becomes unstable to small density fluctuations. An equivalent 
approach is the Vlasov formalism, a semi-classical limit of the description of 
the collisionless regime, which has been used in lfT4ll to calculate the dynamical 
spinodal within several RMF models. Contrary to the calculation of the ther- 
modynamical spinodal, the effect of the finite range of the nuclear interaction 
as well as the Coulomb interaction and the presence of electrons is taken into 
account in the determination of the dynamical spinodal. 

In the Vlasov formalism the equilibrium state, characterized by the Fermi 
momenta of neutrons, protons and electrons, ppn, PFp, PFe, is the starting 
point. Charge neutrality requires ppe = PFp- A perturbation of the system 
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is then described by the perturbed mesonic fields, Fj = Fjo + 5Fi , and the 
perturbed distribution functions for the neutrons, protons and electrons 

/i(t,r, p) = foi + Sfi , 6fi = {SiJoi}, 

where /oj = 0{p'jp_ — p^) and the generating function 5(r, p, t) = 
diag (5p, Sn, Se) has been introduced. The time evolution of the distribution 
function fi is described by the Vlasov equation |[T4l 

dfi 

-^ + {fi,hi} = 0, i=p,n,e. 

In the limit of small perturbations the linearized Vlasov equation 

— -^ + {Si,hoi} = 5hi 
at 

is solved. In this equation /iQi and 6hi stand, respectively, for the equilibrium 
single particle Hamiltonian and corresponding perturbation from equilibrium, 
and {v, w} represents the Poisson bracket of two dynamical functions v and w. 
Longitudinal fluctuations are described by the ansatz 

{Si SFj 5 Pi Sh^) = {S^,^{x) 6F^j 5p^,^ 5K,i)e'^'i-'-^'\ 

where q and u are the transferred momentum and energy and x = cos(p • q). 
The dynamical spinodal surface is characterized by a zero frequency uj. 

It is expected that the transition density lies in the metastable region between 
the binodal surface and the dynamical spinodal surface. The binodal surface is 
defined in the p, Yp, T phase space as the surface where the gas and liquid 
phases coexist, and it defines an upper limit for the extension of the pasta phase 
since it also does not take into account Coulomb nor finite size effects. The bin- 
odal surface is calculated imposing Gibbs conditions: for a given temperature, 
the pressure and the proton and neutron chemical potentials are equal in both 
phases ifTTl . The thermodynamical spinodal touches the binodal surface at the 
critical point, which, for a given temperature, occurs for the largest pressure on 
both surfaces and at a density and proton fraction close to the crust-core tran- 
sition density and proton fraction of cold stellar /3-equilibrium matter. This is 
represented in Fig. [2^) where the square and the circle represent respectively 
the crust-core transition from the thermodynamical spinodal and binodal. Both 
the EOS of /3-equilibrium neutrino free matter and matter with trapped neutri- 
nos are represented. We conclude that for neutrino free matter both estimations 
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almost coincide, while for matter with a lepton fraction 1/^=0.4 there is a large 
difference. As stated before, it is possible to get a better estimation of the lower 
limit of the pasta phase extension if, instead of the thermodynamical spinodal, 
the dynamical spinodal is calculated, as in Fig. [2J5). Taking into account finite 
range effects and electrons makes the spinodal region smaller, so even though 
the thermodynamical method gives a good estimation of the transition density 
for cold /^-equilibrium matter, it is a bit too large. This is confirmed by pasta 
calculations f\M as seen in the next section. 
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Figure 2. The transition density for neutrino free /3-equilibrium matter and mat- 
ter with trapped neutrinos and a fraction of leptons Yl = 0.4: a) estimation 
from the spinodal and binodal sections; b) estimation from the dynamical and 
thermodynamical spinodal. 



For densities inside the dynamical spinodal matter is non-homogeneous, and 
for densities outside the binodal surface, matter is homogeneous. Between these 
two surfaces we may find matter in a metastable configuration. The most prob- 
able configuration is the one with the smallest free energy density and requires 
a pasta phase calculation, the topic of the next section. 



3. Cold and Warm Pasta Phase 



In the inner crust of a neutron star (NS), a liquid gas phase transition can give 
rise to the existence of the pasta phase, which is a frustrated system that arises 
from the competition between the strong and the electromagnetic interactions. 
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The pasta phase appears at densities on the order of 0.001 - 0. 1 fm~^ in neutral 
nuclear matter and at a smaller density range in /3-equilibrium stellar matter. 
The basic shapes of these structures were named as droplets (bubbles), rods 
(tubes) and slabs for three, two and one dimensions respectively. The ground- 
state configuration is the one that minimizes the free energy. In what follows 
we use two different prescriptions in order to build the pasta phase at zero and 
finite temperatures: one is based on phase coexistence (CP) and obeys the Gibbs 
conditions, and the other is the Thomas-Fermi (TF) approximation. Within the 
CP method, the pasta structures are built with different geometrical forms in a 
background nucleon gas. This is achieved by calculating from the Gibbs' condi- 
tions the density and the proton fraction of the pasta and of the background gas, 
so that in the whole we have to solve simultaneously various equations. These 
equations are related to the fact that the pressure, proton and neutron chemical 
potentials and temperature are the same in both phases. Two more equations 
are related to the nucleon effective mass in each phase, and an equation that 
balances the amount of protons in each phase has also to be solved: 

fi4) + {l~f){pi') = Y,p, (8) 

where I and II label each of the phases, / is the volume fraction of phase I: 

P^ - P^^ 
where the total baryonic density is 

P = Pp + Pn, (10) 

and Yp is the global proton fraction. 

The correct parametrization of the surface energy, which is temperature, 
proton fraction and geometry dependent, must be used [.17.1 . The following 
functional for the surface tension coefficient, a, is used, 

a = a{x,T = 0)[l-a{T)xT-h{T)T^] , (11) 

where x = 6'^ stands for the square of the asymmetry parameter defined in 
eq.©. The CP method is very easy to implement; however, since it is not self- 
consistent, the results obtained within the method should be taken with care. 
The self-consistent Thomas Fermi results should be compared with other more 
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Figure 3. Free energy for the homogeneous and pasta hke matter. The red 
dashed line defines the pasta free energy. The crossing of the pasta and the 
homogeneous free energy define the extension of the pasta phase. 



realistic, yet more involving methods that do not impose a pre-defined form, 
such as the quantum molecular dynamics calculation in ifTSl or the 3D Skyrme- 
Hartree-Fock method at finite temperature used in 1191 . 

In Fig. [3] we plot the free energy per particle for homogeneous and pasta 
like matter described by the NL3 model within CP, with T=5 MeV and proton 
fractions Y^ = 0.5, 0.3 and for /3-equilibrium matter with trapped neutrinos 
for a lepton fraction Yl = 0.4. This figure illustrates well the effect of the 
pasta phase on the free energy: the formation of a non-homogeneous phase 
is energetically favored. The system in equilibrium chooses the configuration 
with the lowest free energy, so the pasta-like matter defines the ground state of 
the system if its free energy is lower than the free energy of the corresponding 
homogeneous matter. The upper density limit of the pasta phase lies inside 
the binodal surface and decreases when the proton fraction decreases. For /3- 
equilibrium matter this limit defines the NS crust-core transition. 

In Fig. m we compare the estimations for transition density from a non- 
homogeneous phase to a homogeneous phase obtained from the binodal surface 
(Bin) 111], the dynamical spinodal surface (Sp-d) |[T3llT4l the thermodynamic 
spinodal surface (Sp-t) [12] and from the CP and TF methods. The TF results 
lie always between the results obtained from the dynamical spinodal and the 
binodal surfaces. This is a self-consistent method that should satisfy these two 
constraints. The CP calculation may predict results at densities that lie below the 
dynamical spinodal value, namely at high temperatures and symmetric matter. 
For CP to give reasonable results, it is important that a good parametrization of 
the surface energy is used. The thermodynamical-spinodal calculation always 
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Figure 4. Transition density, for several temperatures and proton fractions Yp = 
0.5 (upper plots), 0.3 (lower plots) and /3-equilibrium (eq-/3) at T=0 (left part of 
upper plots), obtained using different methods within the models a) NL3 and b) 
TW.Fig. taken from m. 



gives a quite good estimation even though it does not take into account neither 
the surface nor the Coulomb effects. For /3-equilibrium matter all methods give 
similar results. This is due to the occurrence of the transition density close to 
the critical point where the spinodal and binodal surfaces touch and the pressure 
on these surfaces is maximum. 

In Fig. \5\ we compare the density range for which each pasta configuration 
exists within NL3 and TW. The thick lines stand for the Thomas-Fermi calcu- 
lation, the thin ones for the CP method. Full lines represent TW and dashed 
ones NL3. For symmetric matter the main difference between the models is 
the appearance of the different phases at slightly smaller densities within NL3. 
However, the largest differences occur for Yp = 0.3: NL3 has a much smaller 
rod phase and no slab phase at T-1 and 8 MeV and quite large tube and bub- 
ble phases. The figure also illustrates the power and limitations of the simpler 
method CP: the onset at lower densities of the pasta phase is quite well de- 
scribed, however, it fails to describe the bubble phase and predicts a smaller 
crust-core transition density. 

It is important to point out that Figs. |4] and [5] are slightly different from 
the ones published in 1 16] for the cases Y^ = 0.3 and matter in /3-equilibrium 
because a slightly different paramerization of the surface coefficient was used, 
i.e., in eq.©, the proton fraction was taken as the proton fraction of the denser 
phase in |[T6i and here it was taken as the global proton fraction. 
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Figure 5. Pasta phases: comparison between NL3 (dashed line) and TW 
(full line) and the methods CP (thin red lines) and TF (thick lines) for a) /?- 



equilibrium cold stellar matter; b) Yp 
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0.5, c) Yp = 0.3. Fig. taken from 



The importance of the a particles cannot be neglected. It is the most strongly 
bound system among all light clusters and it certainly plays a role in nuclear 
matter, mainly at low temperatures. The Lagrangian density that describes the 
a particles can be added to Eq. [T]and is given by |[20i : 



Ca = -{iD^ct^aTiiD^a^ 



-(t>*M: 



a'rai 



(12) 



id^ - r„y , where M^ = 4M - Ba, Ba = 28.3 MeV. 



with iDc 

The coupling of the uj meson to the a-particles is included for mimick- 
ing the excluded volume effect and, consequently, the a particles dissolution 
at high densities. The dissolution density obtained using the ansatz F^ = 4r„ 
is in agreement with the dissolution densities obtained within a quantum statis- 
tical approach [20]. More careful studies are necessary in order to determine 
the adequate meson-cluster couplings. When a particles are included, Eq.® 
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(13) 



becomes 

fip'j, + 2pi) + (1 - f){pl' + 2p'J) = Ypp, 

where I and II label each of the phases, / is the volume fraction of phase I given 
in eq.©, where the total baryonic density is now 



P = Pp + Pn + 4:Pa, 

and Yp is the global proton fraction given by 



Y„ 



Pp + 2/3q, 



P 



(14) 



(15) 



The a particle densities are plotted in Fig. [6]for Yp = 0.5 and 0.3 and T = 5 
and 8 MeV. We include the calculation for both homogeneous matter and pasta- 
like matter. This figure gives a hint on the possible effects of a particles in the 
inner crust of a compact star, which is larger for the larger temperatures and 
larger proton fractions. Due to the existence of a non-homogeneous phase, the 
a-particles dissolve at larger densities, although the a fraction may take very 
small values. It is important to include other small clusters (deuteron, tritium 
and helium 3), which, due to their smaller masses will predominate over the a 
particles at the larger temperatures. We expect that the appearance of these light 
clusters will affect heat and transport properties. 
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10 " Id" 
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Figure 6. a particle density for a) T=5 MeV and b) T= 8 MeV and Yp = 0.3 
and 0.5 obtained with NL3 for homogeneous matter (thin lines) and gas phase 
of the pasta-like matter (thick lines). Fig. taken from 1.17.1 . 

The properties of the pasta phase are also discussed in the chapter Nuclear 
pasta in supemovae and neutron stars by G. Watanabe and T. Maruyama, where 
different formalism and numerical methods are used. 
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4. Constraints on the Equations of State 

We have seen that the EOS that describe equally well the properties of nuclear 
saturation matter and the ground-state properties of nuclei, predict quite differ- 
ent spinodal surfaces at high asymmetries and/or temperatures in Fig [T] The 
same occurs for other properties at high densities, such as the incompressibility. 

The asymmetric nuclear matter EOS may be constrained by the properties of 
asymmetric nuclear matter obtained from various analyses of experimental data, 
including isospin diffusion measurements 11211 . giant resonances ||221, isobaric 
analog states [23] or meson production (pions |[24ll . kaons |[25l ) in heavy ion 
collisions (see |26| for an overview). 

Some correlations between finite nuclei properties and bulk matter proper- 
ties have been obtained: a linear relation between the density derivative of the 
neutron matter EOS at 0.1 fm^^ and the neutron skin thickness of ^ospb that has 
been theoretically tested with different Skyrme parameter sets ll27l and relativis- 
tic Hartree models [28|. Another linear correlation was obtained between the 
208p|j skin thickness and the liquid-to-solid transition density in neutron stars 
||2911. In what follows the neutron skin thickness is defined as 



6R = Rn - Rp, (16) 

where the mean square radius is 

J d-irpi{Y) 

Accurate measurements of neutron skin thicknesses, via future parity violating 
experiments [30 1 or by means of existing antiprotonic atoms data 131 , 321 . are 
thus helpful in determining the bulk properties of nuclear systems. 

The other important quantity, the slope of the symmetry energy, is given by 



dp ^'='">' '^"'~ 2 86^ 






(18) 

<5=0 



Using a wide range of effective nuclear models, both phenomenological 
such as the non-relativistic Skyrme forces or RMF models and microscopic 
Brueckner-Hartree Fock with the realistic AVIS potential plus a three-body 
force of the Urbana type [33], we have confirmed the existence of a linear corre- 
lation between the symmetry energy slope L and the crust-core transition den- 
sity obtained from the thermodynamical spinodal, see Fig. |7^). The same mod- 
els predict a linear correlation between the neutron skin thickness of the ^ospj^ 
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Figure 7. Correlations between a) the symmetry energy slope L and the tran- 
sition density obtained from the thermodynamical spinodal; b) the symmetry 
energy slope L and the neutron skin thickness for "^^^Pb; c) the neutron skin 
thickness for "^^^Pb and the the transition density. Fig. taken from ll33l . 

and L and, consequently, a correlation between the neutron skin thickness of 
the 208 pjj ^j^j j-j^g transition density pt, as proposed in 1,29 J . However, in 1.341 
it was shown that no similar correlation exists between the transition pressure 
and the slope L, Fig. [H due to the large dispersion of the predicted transition 
pressure obtained when independent models are considered. This means that an 
experimental determination of L alone is not enough for a good estimation of 
the crust mass and moment of inertia of a compact star, since it is the transition 
pressure that allows a prediction of the EOS from the observation of glitches 

Eoi. 



5. Diffusion Coefficients and the Pasta Phase 

The neutrino signals detected by astronomers can be used as a constraint to infer 
protoneutron star (PNS) composition. For the same purpose, theoretical studies 
involving different possible equations of state obtained for all sorts of matter 
composition have to be done because the temporal evolution of the PNS in the 
so-called Kelvin-Helmholtz epoch, during which the remnant compact object 
changes from a hot and lepton-rich PNS to a cold and deleptonized neutron star 
depends on two key ingredients: the equation of state (EOS) and its associ- 
ated neutrino opacity |[35l . All contributions from neutrino opacities are related 
with the diffusion coefficients and can to be used as input to the solution of the 
transport equations in the equilibrium diffusion approximation to simulate the 



The inner crust and its structure 



15 







' 


(J 


' 






0.6 


oo* 






1 

> 


0.5 
0.4 
0.3 


- • • 


0» 

o 
o 


O - 

.0 ° 

• o •• 




CL 


0.2 

0.1 




o 




o 








40 


60 


80 100 120 








L 


(MeV) 





Figure 8. Transition pressure obtained from the thermodynamical spinodal as a 
function of the symmetry energy slope L for several Skyrme models (full dots), 
RMF models (empty dots) and a Brueckner-Hartree-Fock calculation |[34ll . 



Kelvin-Helmholtz phase of the protoneutron stars. 

The diffusion coefficients that are related to the neutrino opacities are calcu- 
lated in such a way that the formation of nuclear pasta at low densities is taken 
into account. The diffusion coefficients are given by |[35l 



D, 



de, etXuieu)fu{e,){l - /,(e,)), k = 2, 3,4, 



(19) 



where Xui^i^) is the total mean free path of neutrinos and fui^i^) is the Fermi- 
Dirac distribution. 

To calculate neutrino opacities and mean free paths we consider neutral cur- 
rent scattering reactions Uf, + n ^ Ve + n and v^, + p ^ Ve + P and charged 
current absorption reactions u^ + n ^ e^ + p and Ve + p ^- e~^ + n. Ba- 
sically, it consists in calculating the cross sections an,p for neutrino-nucleon 
scattering reactions and the cross section aa for neutrino absorption reactions 
for both nondegenerate and degenerate thermodynamic limits as done in f36\ . 
The thermodynamic regions of intermediate degeneracy are handled by a simple 
interpolation. The mean free path is given by X^i^u) = — ^ , ^^ , ^ , where 
p = Pp + pn and they are related to the diffusion coefficients by 



A 



/o°° d'^v e^/i.(er.)(l - fy{(^y)) 



(20) 



The diffusion coefficients D2, -D3 and D^ are strongly dependent on the 
EOS and are functions of three thermodynamic variables: p, T and Yi. We 
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start by fixing p, T and Yi from the EOS to calculate the cross sections Up, 
cjn and (Ta as function of the neutrino energy and then we integrate in neutrino 
energy. The numerical procedure used to calculate the diffusion coefficients to 
homogeneous and inhomogeneous matter is the same, except for the nucleon ef- 
fective mass. The pasta structure is obtained by the coexistence phases method, 
which is based on the enforcement of the Gibbs conditions. Hence, all other 
thermodynamic variables (chemical potentials of all particle species, tempera- 
ture, pressure and lepton fraction), necessary to calculate the neutrino mean free 
path, and consequently the diffusion coefficients of the pasta phase are equal in 
both phases. The nucleon effective mass, on the other hand, is not the same. To 
calculate Xy{ey) we need the nucleon effective mass M*. In the pasta phase, 
two different phases coexist (phase I and phase II). In our calculation, we have 
used M* = fM<^^ + (1 - f)M<^^^ for the pasta phase. 

Our results for the diffusion coefficients as a function of the baryon density 
at temperature T = 5 MeV and lepton fraction Yl = 0.4 (includes electrons 
and trapped neutrinos) is shown in Fig. |9l from where we can see that only three 
structures are found inside the pasta phase for the present model: droplets, rods 
and slabs. While the diffusion coefficients obtained with homogeneous matter is 
always smooth and continuous, a common trend of all the diffusion coefficients 
obtained with the pasta phase is a kink at very low densities in between 0.01 and 
0.015 fm^^ due to the fact that the effective nucleon mass becomes greater than 
the corresponding chemical potential. 

In obtaining the diffusion coefficients, the EOS was calculated as a grid 
where temperature ranges are in between and 50 MeV and densities vary from 
0.005 to 0.5 fm~^. We have calculated the diffusion coefficients only for bary- 
onic densities above 0.005 fm~^ because the integrals of type (fT9l) are very 
difficult to converge at lower densities. In all cases the diffusion coefficients ob- 
tained with homogeneous matter join the curves obtained with the pasta phase 
at densities higher than the ones shown in Fig |9] For D2 calculated at T=5 
MeV and Yl = 0.4, for instance, they cross each other at p = 0.12 fm~^. Our 
codes interrupt the calculation once homogeneous matter becomes the ground 
state configuration. This means that there is always a gap in the diffusion co- 
efficients when the transport equations are calculated with the inclusion of the 
pasta phase. 

While the diffusion coefficients obtained with homogeneous matter are al- 
ways smooth and continuous, a common trend of all the diffusion coefficients 
obtained with the pasta phase is a kink at very low densities in between 0.01 
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Figure 9. Diffusion coefficients D^, D^ and D^^ as function of baryon density 
at T = 5 MeV and Yl = 0.4 for iiomogeneous matter and pasta phase. Fig. 
taken from [37]. 
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and 0.015 fm"-^. The interpolation procedure we use depends on the quantities 
Vi = if^i ~ M*)/T,i = p,n. Whenever either 7]p or r/„ inverts its sign, these 
kinks appear, i.e., they are the result of the effective nucleon mass being greater 
than the corresponding chemical potential. Moreover, the pasta phase diffusion 
coefficients are always lower than the corresponding coefficients obtained with 
homogeneous matter. 

Our results show that the mean free paths are significantly altered by the 
presence of nuclear pasta in stellar matter when compared with the results ob- 
tained with homogeneous matter These differences in neutrino opacities will 
have consequences in the calculation of the Kelvin-Helmholtz phase of pro- 
toneutron stars 1.37.1. 



6. Quark Stars Subject to Strong Magnetic Fields 

Neutron stars with very strong magnetic fields of the order of 10^^ — 10^^ G 
are known as magnetars and they are believed to be the sources of the intense 
gamma and X rays detected in 1979 |[38l[39l . The hypothesis that some neutron 
stars are constituted by unbound quark matter cannot be completely ruled out 
since the Bodmer-Witten conjecture |T| cannot be tested in earthly experiments. 
This conjecture implies that the true ground state of all matter is (unbound) 
quark matter because theoretical predictions show that its energy per baryon at 
zero pressure is lower than ^®Fe binding energy. According to 131 the quarks 
are bound by the strong force rather than the gravitational force that binds other 
stars. At the surface of the star, where the quark density drops abruptly to zero, 
the electrons extend into a layer of thickness of the order of 10^ fm above the 
surface, where there is a super-strong electric field that ties the electrons to the 
star If such a star would be covered by a crust of ordinary nuclear matter (and 
neutralizing electrons), it would be blown away as it forms Ii40il or would be 
destroyed by thermal effects [41] . Thus, one expects quark stars to be bare, in 
the sense that the surface is this thin layer of electrons ll42l . 

We investigate quark matter described by the Nambu-Jona-Lasinio H3]| 
model exposed to strong magnetic fields. In the case of pure quark matter, as 
predicted by the QCD phase transition possibly taking place in heavy ion colli- 
sions, the ultra-strong magnetic field results from the superposition of external 
and internal fields, the former generated by the alignment of charged particles 
that are spinning very rapidly. We use an external field to mimic the real situa- 
tion. 
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In order to consider (three flavor) quark matter subject to strong magnetic 
fields we introduce the following Lagrangian density 

C = Cf- ^F^uF^"" (21) 

where the quark sector is described by the su(3) version of the Nambu-Jona- 
Lasinio model which includes scalar-pseudoscalar and the t'Hooft six fermion 
interaction that models the axial U{1)a symmetry breaking: 

Cf = V-/ [7a. {id'' - QfA'^) - rhc] Ipf + Csym + Cdet , (22) 

where Csym and Cdet are given by: 

8 

Csym = GY, [(V^/AaV'/)' + (V^/nsAaV-/)'] , (23) 

a=0 

Cdet = -K {detj [VS/(1 + 75)^^/] + det/ [ijf{l - 75)^'/] } , (24) 

where V/ = {u, d, s)^ represents a quark field with three flavors, rhc = 
diag f{'mu,mii,'ms) is the corresponding (current) mass matrix while qf rep- 
resents the quark electric charge, Aq = a/2/3/ where / is the unit matrix in the 
three flavor space, and < Aa < 8 denote the Gell-Mann matrices. Here, we 
consider niu = rud 7^ rus- A^ and F^jy = df^A,^ — d^A^ are used to account 
for the external magnetic field. Since we are interested in a static and con- 
stant magnetic field in the z direction, A^ = 6^2XiB. Whenever /3-equilibrium 
neutrino-free matter is considered, the leptonic sector is given by eq.®, where 
neutrinos are not taken into account. Besides the su(3) version of the model, we 
also present two flavor su(2) NJL model results for comparison. 

In Fig. [TOk ) one can see that the inclusion of the magnetic field makes 
matter more and more bound in both versions of the model. For the present set 
of parameters, the energy per baryon E/A of magnetized quark matter becomes 
more bound than nuclear matter made of iron nuclei, ^Isep^ ~ 930 MeV for B 
around 2 x lO^^ G. 

The EOS for stellar matter under a strong magnetic field is obtained with 
a density dependent frozen magnetic field which is set equal to 10^^ G at the 
surface and does not exceed 6 x 10^® G in the center of the star. Approximate 
values of the mass and radius of the stars obtained from the integration of the 
TOV equations [i45il are displayed in Fig. [TOb). from where it is seen that the 
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Figure 10. a) Energy per nucleon as a function of density for B = 0, 10^^ 
and 2 x 10^^ G within NJL su{2) and NJL sn(3); b)Mass-radius curves for the 
families of stars within NJL su{2) and NJL su{3). For B < 10^^ the curves 
coincide with the B = results. Fig. taken from 1,44,1 . 



gravitational mass increases with the increase of the magnetic field for an in- 
tensity larger than '--^ 5 x 10^^ G for the su(3) version and 10^® G for the su(2) 
NJL. Another important effect of the field on the properties of the stars is the 
increase of the radius of the star. The largest radius may be as high as 9.5 Km 
for the su{3) NJL. In general, the maximum mass star configurations for the 
su{2) version of the NJL model are smaller with smaller radius, ~ 7 Km. 

Next we discuss the effect of the magnetic field and temperature on quark 
matter as possibly formed in heavy-ion reactions. We display the free energy per 
particle in terms of the magnetic field for symmetric matter in Fig. fTTk ) and in 
terms of the density for asymmetric matter in FigJTTb). As discussed before, the 
effect of the magnetic field is stronger for smaller densities and temperatures. 
At a given density the main effect of temperature is to decrease the free energy 
per particle. The effect of the magnetic field is clearly seen in Fig. ITTh). It is 
stronger for the smaller temperatures and, for a strong enough field it gives rise 
to a decrease of the free energy (above B = A x 10^^ G for p = po), due to 
a reduction of the number of Landau levels. However, for even stronger fields, 
the free energy increases due to an increase of the effective quark masses with 
B. This explains why in FigJTTb') the free energy at T=10 MeV below p = 0.5 
fm^'^ is smaller for S = 2 x 10^^ G and larger for B = 10^° G. For larger 
temperatures, the reduction observed at intermediate densities washes out and it 
is only observed an increase of the free energy for fields above 5 x 10^^ G. The 
increase of the free energy for very large fields is mainly due to a reduction of 
the entropy and an increase of the effective mass for B > 10^^ G. 
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Figure 1 1. Free energy per particle for NJL a) su(3) at T = 10, 100, 150 MeV 
as a function of B for symmetric quark matter and p = po and 3po ; b) su(2) 
versus the density for asymmetric matter pd = l.2pu, T = 10, 100, 150 MeV 
and B = 10^*^, 2 x lO^^, 10^° G. Fig. taken from |i46J 

To extend the present calculation in order to consider matter in protoneu- 
tron stars, the enforcement of charge neutrality and /3-equilibrium is required. 
The ideal calculation would consider a fixed entropy resulting in temperatures 
lower than 50 MeV in the interior of the stars. This extension is simple and 
straightforward. 



7. Conclusions and open questions 

In this chapter we have focused on the study of possible internal structures of 
the neutron stars crust. The existence of the pasta phase is related to the in- 
stabilities and possible coexisting phases in the nuclear matter, features well 
described by the spinodal and binodal sections. Important remaining questions 
are: at which temperatures do the pasta structures dissolve? Up to which tem- 
peratures is the ansatz of a Wigner Seitz cell valid? We have tried to validate 
a simple coexistence phase approach to obtain the pasta phase and its structure 
by using the results of possible surface energy parametrizations and compar- 
ing the final results with the more realistic and self-consistent Thomas-Fermi 
approximation, which is very time consuming. It remains to be checked how 
the Thomas-Fermi calculation of /3-equilibrium pasta phase compares with the 
results obtained within Hartree-Fock and Hartree-Fock-Bogoliubov, where only 
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spherical clusters are considered. More calculations in this direction are still 
required. 

So far, we have included a particles as part of the pasta phase. Will light 
clusters other than a particles contribute in a significative manner to the ex- 
tension of the pasta phase and its transition density to homogeneous matter? 
Moreover, how will the presence of light clusters affect the properties of the 
pasta phase? Are the values for the dissolution density of the different light 
clusters in nuclear matter realistic BTl ? 

We have also seen that the inclusion of the pasta phase affects the neutrino 
opacities through its diffusion coefficients. It remains to be checked how the 
internal structure of the pasta phase, which is model and approximation depen- 
dent, influences the diffusion coefficients through the different neutrino mean 
free paths. A more complete calculation will show us how the existence of the 
pasta structures affect the transport properties of the crust. 

Our last section was devoted to a subject other than the crust of neutron 
stars. The measurement of a pulsar of inferred mass 1.97 ± 0.04Mq known as 
PSR J1614-2230 HSl . opened new questions about its possible constituents and 
related equation of state. It is well known that there is a maximum mass that can 
be supported against collapse and its theoretical value is very model dependent. 
According to P9l . the maximum mass of a quark star lies between 2 and 2.7 
Mq. Hence, it is claimed that PSR J 16 14-2230 can be a quark star. Do quark 
stars really exist? How are they affected by strong magnetic fields? We have 
seen that our results for the NIL model cannot reproduce such high masses, not 
even with very strong magnetic fields. 
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